Are external perturbations responsible for chaotic motion in galaxies? 
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Abstract 

We study the nature of motion in a logarithmic galactic dynamical model, with an additional external perturbation. 
Two different cases are investigated. In the first case the external perturbation is fixed, while in the second case it 
is varying with the time. Numerical experiments suggest, that responsible for the chaotic phenomena is the external 
perturbation, combined with the dense nucleus. Linear relationships are found to exist, between the critical value 
of the angular momentum and the dynamical parameters of the galactic system that is, the strength of the external 
perturbation, the flattening parameter and the radius of the nucleus. Moreover, the extent of the chaotic regions in 
the phase plane, increases linearly as the strength of the external perturbation and the flattening parameter increases. 
On the contrary, we observe that the percentage covered by chaotic orbits in the phase plane, decreases linearly, as the 
scale length of the nucleus increases, becoming less dense. Theoretical arguments are used to support and explain the 
numerically obtained outcomes. A comparison of the present outcomes with earlier results is also presented. 
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1. Introduction 

Most galaxies in the Universe, are gravitationally bound 
to a number of other galaxies. These form a fractal-like 
hierarchy of clustered structures, with the smallest such 
association being termed groups. A group of galaxies is 
the most common type of galactic cluster and these for- 
mations contain the majority of the galaxies. In other 
words, galaxies interact with each other, most common 
in pairs or as triple systems. The Milky Way-Magellanic 
Clouds system, is a very good example of galactic interac- 
tion [1,8,10,16,20,24,27,28]. Another interesting system of 
interacting galaxies, is the Andromeda M3I galaxy, with 
its two small companion galaxies M32 and NGC 205. This 
interesting triple system was investigated, in an earlier pa- 
per, using a self consistent computer simulation code, with 
interesting results [26]. Moreover, the formation of spiral 
structure in a galaxy, as a result of the gravitational per- 
turbation caused by a permanent companion, was studied 
in a previous work [25]. 

In this article, we shall use the gravitational potential 



V (r,z) = ^-ln [r 



(1) 



in order to describe the motion in our galactic dynamical 
system. This potential is important for galactic dynamics 



and represents an elliptical galaxy, with a nucleus of ra- 
dius c, which displays a flat rotation curve at large radii 
[2]. Here r,z are the usual cylindrical coordinates. The 
flattening parameter I < a < 2, defines the axial ratio of 
the equipotential ellipsoids. Logarithmic potentials, have 
been frequently used by many researchers, over the last 
decades, in order to model galactic motion [6,12,19,21,22]. 
The parameter 0.01 < c < 0.25 is the radius of the nu- 
cleus, while vq is a parameter used for the consistency of 
the galactic units. In order to keep things simple, we con- 
sider that our galaxy is subject to an external perturbation 
V\ = — Ar, caused by a companion galaxy, where A > is 
the strength of the external perturbation. Thus the total 
potential is 



V(r,z) = V + V 1 = ^-\n[ 



r 2 + az 2 + c 2 ] - Ar. (2) 
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We use a system of galactic units, where the unit of 
length is lkpc, the unit of time is 0.97746 x 10 8 yr and 
the unit of mass is 2.325 x 1O 7 M . The velocity and the 
angular velocity units are 10km/s and lOkm/s/kpc respec- 
tively, while G is equal to unity. The energy unit (per unit 
mass) is 100(km/s) 2 . In the above units, we use the value 
vq = \bkm/s, while a, A and c are treated as parame- 
ters. We consider only bounded motion, therefore we take 
A < 21. 

As the total potential V(r, z) is axially symmetric and 
the L z component of the angular momentum is conserved, 



Published in Chaos, Solitons & Fractals 



March 2, 2013 



we use the effective potential 

Veff(r,z) = ^ + V(r,z), (3) 

in order to study the character of motion in the meridian 
r — z plane. The equations of motion are 

„ = dV ef f(r,z) 
dr 



where the dot indicates derivative with respect to time. 
The corresponding Hamiltonian is written as 

H=±(p 2 r +p 2 z )+V eff (r,z)=E, (5) 

where p r andp z are the momenta, per unit mass, conjugate 
to r and z respectively, while E is the numerical value of 
the test particle's energy, which is conserved. Equation 
(5) is an integral of motion, which indicates that the total 
energy of the test particle (star), is conserved. 

The outcomes of the present research are based on the 
numerical integration of the equation of motion (4). We 
use a Bulirsh-Stoer integration routine in Fortran 95, with 
double precision in all subroutines. The accuracy of our 
results was checked by the constancy of the energy inte- 
gral (5), which was conserved up to the twelfth significant 
figure. 

Our aim is to connect the external perturbation caused 
by a companion galaxy, with the character of motion (reg- 
ular or chaotic). Moreover, we try to find relationships, 
in order to connect the dynamical parameters of the sys- 
tem, that is the strength of the external perturbation, the 
flattening parameter and the radius of the nucleus, with 
the chaotic percentage and the critical value of the angular 
momentum L zc . Before doing this, it would be interesting 
to try to relate, the chaotic orbits or the chaotic motion, 
in general, with actual astronomical issues. 

One of the main reasons for the transition from or- 
dered to chaotic stellar motion, is the presence of a mas- 
sive object in the central regions of galaxies. Stars reaching 
the center on highly eccentric radial orbits, are scattered 
out of the galactic plane, displaying chaotic motion [23]. 
Furthermore, a central mass concentration can strongly 
perturb the stellar orbits in elliptical galaxies, which be- 
come chaotic [18]. Observational indications suggesting, 
the presence of strong central mass concentration, with 
a very sharply raising rotation curve. A second reason, 
for chaotic behavior in galactic motion, is the presence 
of strong external perturbations. As the perturbation in- 
creases, it destroys the stability of orbits and increases 
the amount of the stochasticity of the dynamical system. 
Resonances are also responsible for chaotic motion [9,11]. 
Here we must emphasize, that resonances are not only re- 
sponsible for the presence of chaos in galaxies, but also for 
the chaotic motion in the solar system [14,29]. The reader 
can find interesting information on the chaotic motion in 



galaxies and its connection with observations, in the work 
of Grosb0l [13]. 

The layout of this article is organized as follows. In 
Section 2, we study the structure of the phase plane and 
derive numerically relationships between the percentage 
of the area A% covered by the chaotic orbits in the phase 
plane and the parameters A, a and c. In the same Sec- 
tion, we present relationships between the critical value of 
the angular momentum and the same basic parameters of 
the dynamical system. In Section 3, these relationships 
are also reproduced and explained, using theoretical argu- 
ments. In Section 4, we follow the evolution of the orbits, 
as the flattening parameter, or the strength of the exter- 
nal perturbation changes with time. We close with Section 
5, where a discussion and the conclusions of the present 
research are given. 

2. Structure of the dynamical system 

Figure la-d shows the r — p r ,z — 0,p z > Poincare 
phase plane, for the Hamiltonian (5) obtained by numer- 
ical integration of the equations of motion (4). Figure la 
shows the r — p r phase plane when: A = 0, a — 1, L z = 10 
and E = 467. Here we have the case of a spherical un- 
perturbed galaxy and therefore, as we expected, the entire 
phase plane is occupied by invariant curves produced by 
regular orbits of the 1:1 resonance. Figure lb shows the 
r — p r phase plane when: A = 21, a = 1,L Z = 10 and 
E = 300. In this case, the spherical galaxy is affected by 
the largest permissible value of the external perturbation 
A. Here we observe, that the majority of orbits are reg- 
ular orbits. The external perturbation has as a result a 
considerable chaotic layer, which is confined in the outer 
parts of the phase plane. Secondary resonances, if any, are 
negligible. Figure lc shows the r — p r phase plane when: 
A = 0,a = 1.9, L z = 10 and E = 467. Here we have the 
case of the unperturbed, non spherical galaxy. As one can 
see, again the majority of orbits are regular, while a very 
small chaotic layer in the outer parts of the phase plane 
is present. Figure Id shows the r — p r phase plane when: 
A = 21, a = 1.9, L z = 10 and E = 300. One observes a 
large chaotic sea, while there are also areas of regular mo- 
tion. In addition to the above, one can also observe smaller 
islands of invariant curves, embedded in the chaotic sea, 
which are produced by secondary resonances. The values 
of all other parameters are: vo = 15 and c = 0.25. The 
values of the energy E, were chosen so that in all phase 
planes r max ~ 8. 

The main conclusion from the above analysis, is that 
spherical galaxies display less chaos, than flat ones, when 
exposed to external perturbations. In other words, we see 
that chaos decreases when the symmetry of the galaxy 
increases. 

Figure 2a-h shows eight representative orbits of the dy- 
namical system. Figure 2a shows a regular orbit when: 
A = 0, a = 1 and E = 467. Initial conditions are: r = 
1.5, z = 0,p r o = 0, while the value of p z o is found from 
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Figure 1: (a-d): The r — p r , z = 0,p z > Poincarc phase planes, for the Hamiltonian (5) when (a-upper left): A = 0, a = 1, E = 467, (b-uppcr 
right): A = 21, a = 1, E = 300, (c-down left): A = 0, a = 1.9, E = 467 and (d-down right): A = 21, a = 1.9, E = 300. The values of the other 
parameters arc: vq = 15, c = 0.25 and L z = 10. 



the energy integral (5), for all orbits. Figure 2b shows a 
regular orbit when: A = 21, a = 1 and E = 300. Initial 
conditions are: r = 2.8,2:0 = 0,p r o = 12- In Figure 2c 
a quasi periodic orbit, characteristic of the 2:3 resonance, 
is shown. Here: A = 21, a = 1.9 and E = 300. Ini- 
tial conditions are: r — 6.3, z — 0,p r o — 0. Figure 2d 
shows a quasi periodic orbit when: A = 0,a = 1.9 and 
E = 467. Initial conditions are: ro = 6.2, zq — 0,p r o = 0. 
This orbit is characteristic of the 4:3 resonance. In Figure 
2e a quasi periodic orbit, characteristic of the 4:5 reso- 
nance, is presented. Here: A = 21, a = 1.9 and E = 300 
while, initial conditions are: r = 0.6, z = 0,p r o = 0. 
In Figure 2f we see a characteristic orbit of the 6:7 reso- 
nance. Here: A = 21, a = 1.9 and E = 300 while, ini- 
tial conditions are: r = 3.5, z — 0,p r o — 0. Figure 2g 
shows a complicated orbit produced by the 8:7 resonance, 
when: A = 21, a = 1 and E = 300. Initial conditions are: 
r = 6.23, z = 0,p r o = 0. A chaotic orbit is given in 
Figure 2h, when: A = 21, a = 1.9 and E = 300. Initial 
conditions are: ro = 0.42, zq — 0,p r o = 0. All orbits were 
calculated for a time period of 100 time units. The values 
of all other parameters are: vq = 15, c = 0.25 and L z = 10. 



It is evident, that responsible for the chaotic regions 
is the external perturbation. This can be justified looking 
at Figure 3a, which displays the percentage of the area 
A% covered by chaotic orbits, in the r — p r , z = 0,p z > 
phase plane as a function of the external perturbation A. 
The values of all the other parameters are: E — 300, a = 
1.9, c = 0.25 and L z = 10. We see that the chaotic area 
A% increases linearly with A. Furthermore, the numerical 
calculations suggest that, for the above given values of the 
dynamical parameters, the chaotic regions are negligible 
when A — > 0. Let us now come to see how the percentage 
of the chaotic regions in the r—p r phase plane, is connected 
with the flattening parameter a. The results are given in 
Figure 4a, when 1 < a < 2. The values of all the other 
parameters are: E — 300, A = 21, c = 0.25 and L z = 10. 
One observes, that the chaotic area A% increases again 
linearly with a. On the contrary, as we observe in Figure 
5a, the percentage of the chaotic regions in the r — p r 
phase plane, decreases this time, linearly as the radius of 
the nucleus increases. The values of all other parameters 
are: E = 300, a = 1.9, A = 21 and L z = 10. Therefore, 
we conclude that galaxies with less dense nucleus, display 
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Figure 3: (a-b): (a-left): Plot of the percentage of the area A% in the phase plane covered by chaotic orbits vs. A, when L z = 10, (b-right): 
Relationship between the critical value of the angular momentum L zc and the external perturbation A. The values of other parameters are: 
a = 1.9, c = 0.25 and E = 300. 
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Figure 4: (a-b): (a-left): Plot of the percentage of the area A% in the phase plane covered by chaotic orbits vs. a, when L z = 10, (b-right): 
Relationship between the critical value of the angular momentum L zc and the flattening parameter a. The values of other parameters are: 
A = 21, c = 0.25 and E = 300. 
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Figure 5: (a-b): (a-left): Plot of the percentage of the area A% in the phase plane covered by chaotic orbits vs. c, when L z = 10, (b-right): 
Relationship between the critical value of the angular momentum L zc and the radius of the nucleus c. The values of other parameters are: 
a = 1.9, A = 21 and E = 300. 



less chaotic motion. An explanation of these numerically 
found relationships of the dynamical system, will be given 
in the next Section. 

Before closing this Section, we would like to present 
relationships connecting the critical value of the angular 



momentum L zc (that is the maximum value of the angular 
momentum L z , for which the test particle displays chaotic 
motion) and the basic parameters of the dynamical sys- 
tem, that is the strength of the external perturbation A, 
the flattening parameter a and the radius of the nucleus 
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c. Our numerical experiments show that the relationship 
between L zc and A is linear. This linear relationship is 
shown in Figure 3b. In Figure 4b, we observe that the re- 
lationship connecting the flattening parameter a with the 
critical value of the angular momentum L zc , ia also linear. 
Orbits starting in the upper part of the a — L zc plane arc 
regular, while orbits starting in the lower part of the same 
plane, including the line display chaotic motion. Figure 
5b shows the relationship between the the critical value 
of the angular momentum L zc and the radius of the nu- 
cleus c. Here one can observe a straight line fitting all the 
numerically found data points. Orbits with values of L zc 
and c on the left part of the c — L zc plane, including the 
line, are chaotic, while orbits with values of the parameters 
on the right part of the same plane are regular. What is 
more interesting, is that all these linear relationships can 
be also derived, using some elementary semi-theoretical 
arguments. 

3. Semi-theoretical arguments 

In this Section, wc shall present some semi-theoretical 
arguments together with elementary numerical calculations, 
in order to explain the numerically obtained relationships 
given in Figs. 3, 4 and 5. 

Our analysis takes place near the nucleus because, there, 
all internal forces acting on the star take their maximum 
values. Furthermore, the tangential velocity also takes 
its maximum value there. Note that the external force is 
constant, for all parts of the galaxy taking the value 

F ext = ~ = A. (6) 

As the test particle approaches the nucleus there is a 
change in its momentum given by 

mA Pr =< F mt + F ext > At =< F int + A > At, (7) 

where m is the mass of the test particle, F^ n t and F ex ± 
is the internal and external average force acting on the 
particle, while At is the duration of the encounter. We 
assume that the star displays chaotic motion, when the 
total change in the momentum after n > 1 encounters is 
of order of 

L zc 



<r > 
Thus we have 



, r « 1. 



(8) 



m 



a p™ «< F mt + \>j2 At *- 

i=l i=l 
Setting 

n 

, Ap ri = mv^ 



(9) 



mL, 



m 



i=i 



<r > 



1=1 



m 



(10) 



in equation (9) we find 

F 7 r- 



< Fint + A > T c 



(11) 



< r > 

Rearranging, we write equation (11) in the form 

< L zc >w L zc0 +mX, (12) 

where jjL\ = T c tq,L zc q = Fi nt T c rQ. As the internal force, 
at a given point near the nucleus (r = r <C c, z = z <C c), 
is fixed, equation (12) gives the linear relationship between 
A and L zc , which explains the numerical results shown in 
Figure 3b. 

The linear relationship between the radius of the nu- 
cleus c and L zc can explained similarly, using again semi- 
theoretical arguments. In fact, we use essentially simi- 
lar arguments to those used in Caranicolas & Innanen [3], 
Caranicolas & Papadopoulos [5], Caranicoals & Zotos [7] 
and Zotos [30]. When the star approaches the dense nu- 
cleus, its momentum in the z direction changes according 
to the equation 



mAp z =< F z > At, 



(13) 



where m is the mass of the star, < F z > is the total average 
force acting in the z direction, while At is the duration of 
the encounter. It was assumed, that the star's deflection 
into higher values of z, proceeds in each time cumulatively, 
a little more, with each successive pass by the nucleus and 
not with a single tragic encounter. It is also assumed, that 
the star is scattered off the galactic plane, after n(n > 1) 
encounters, when the total change in the momentum in the 
z direction, is of order of the tangential velocity v^, of the 
star near the nucleus, at an average distance r =< r >. 
Thus we have 



m 



^Ap zl ^<F Z >J2 A U- 



(14) 



i=i 



If we set 



rn 



<r > 

]TAt, = T e , 

i=i 

m = I, (15) 
and combine equations (14) and (15), we find 



i=i 



< r > 



t< F z > T e . 



(16) 



The force acting in the z direction for a star of unit mass 
(m = 1) is 



dVo —cxvqZ 
dz r 2 + az 2 + c 2 ' 



(17) 
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Setting z = c,r =< ro >~3> z (remember that the star 
before scattering is very close to the galactic plane) and 
keeping only linear terms in z and c, Eq. (17) becomes 



<F Z > r 



<r > 



2 ' 



(18) 



Inserting the value from Eq. (18) into relation (16) we 
obtain 



<r > 



-ap,2C, 



(19) 



where jj,2 = T e / < ro >. Here we must note that, equation 
(19) cannot be considered as an exact representation of 
the relation, between the involved quantities. It rather 
can be seen, as an indication of the relation that needs 
to be completed with additional terms. Those terms, can 
be derived through numerical calculations, providing the 
necessary information. Actually, numerical calculations, 
suggest that expression (19), needs to be supplemented 
with an additional constant term giving the value of L zc 
in the case when c — > 0. Calling this term L zc0 , expression 
(19) takes the form 



< Lzc >~ L zc0 - a/x 2 c. 



(20) 



Relation (20) explains the linear relationship between 
the radius of the nucleus c and L zc , shown in Fig. 5b. 
Equation (20) also contains the flattening parameter a. 
As one can see, for a given value of c, there is a linear 
dependence between a and L zc . This explains the numer- 
ically obtained linear dependence shown in Fig. 4b, where 
c = 0.25, A = 21 and E = 300. 

As the scattering occurs near the nucleus, it must be 
^ < c and zq <C c. The particular values of ro and zq 
are irrelevant. As we can see from Eq. (17), \F Z \ increases 
linearly with the flattening parameter a. As a consequence 
of this linear increase of the \F Z \, we see that the chaotic 
region in the r — p r phase plane, increase to high values, 
up to about A = 82%, when a reaches 2 (see Fig. 4a). 

4. Evolution of orbits in the time-dependent model 

In this Section, we shall study the evolution of orbits as 
the parameters a and A change linearly with time following 
the equations 



a(t) = a in - kit, 
A(i) = A m + kit, 



(21) 



where oti n ,\i n are the initial values of a and A, while k\ 
and &2 are parameters. 

Figure 6a-b shows the evolution of an orbit, as a changes 
with time, following the first of equations (21). The initial 
conditions arc: ro = 4, zq = p r o = 0. The initial value 
of energy is E = 300, A — 21, c — 0.25, while the value 
of p z o is found from the energy integral (5) in all cases. 
The initial value of a is 1.9, while k\ is equal to 0.01. The 
values of all other parameters are as in Fig. Id. Figure 6a 



shows the orbit for the first 90 time units, while figure 6b 
shows the orbit for the rest 110 time units. One observes, 
that the orbit starts as chaotic and tends to be a regular 
orbit as a approaches unity. At t = 90 when a = 1, the 
time evolution stops and the system is now spherical, but 
always subject to a large external perturbation (A = 21). 
The energy value was settled to the value E = 276.48. The 
evolution of the orbit in the spherically symmetric system 
is shown in Figure 6b. As expected, the orbit is now reg- 
ular and remains regular. Of course, the value of energy 
E must correspond to r < 6, for which regular orbits do 
exist. In Figure 6c one observes the complete evolution of 
the orbit for a time period of 200 time units. Figure 6d 
shows the evolution of the Lyapunov Characteristic Expo- 
nent (L.C.E) [17], for a time period of 10 5 time units. For 
the first 90 time units, the profile of the L.C.E indicates 
chaotic motion, but for the rest time interval, the L.C.E 
corresponds to regular motion. Here we must point out, 
that the time scale with which the orbits become regular 
depends on the initial values of the quantities A and c. 
Numerical calculations not given here suggest that, as the 
system evolves from a flat system to a spherical one, the 
majority of chaotic orbits become regular. 

Figure 7a shows the evolution of an orbit, as the ex- 
ternal perturbation A changes with time following the sec- 
ond of equations (21). The initial conditions arc: r = 
2.0, zo = p r o = 0, a = 1.9, c = 0.25, while the initial value 
of energy is E = 300. The value of angular momentum 
is L z — 5. The initial value of A is 0, while hi is equal 
to 0.1. Figure 7a shows the orbit for the first 210 time 
units, while Figure 7b shows the orbit for the rest 90 time 
units. One observes, that the orbit starts as a regular 
orbit and gradually becomes chaotic, as the value of the 
external perturbation A increases. After 210 time units, 
the value of A becomes 21. At this point, the evolution 
of our dynamical system stops and the orbit runs for 90 
more time units, with its new value of energy, which is 
now E = 275.38. In Figure 7c we present the evolution 
of the total velocity v(t) = y/p* + p 2 z of the test parti- 
cle, as a function of time for the above orbit and for time 
interval of 500 time units. One observes that, at about 
t = 200 time units, the velocity profile changes, while at 
the same time the velocity increases. This indicates that, 
in practice, the orbit becomes chaotic after the external 
perturbation has reached the value A = 20. Here one must 
notice, that the above results are in agreement with ob- 
servational data, where an increase of the stellar velocity 
is expected, in regions with significant chaos. Moreover, 
observations show, that in chaotic regions one expects to 
get an asymmetric velocity profile [13]. As one can see in 
Figure 7c, our velocity profile becomes asymmetric, when 
the motion changes from regular to chaotic. In order to 
double check our results, we computed the maximal Lya- 
punov Characteristic Exponent (L.C.E), for a time period 
of 10 5 time units, which is shown in Figure 7d. The L.C.E 
indicates regular motion for the first 210 time units and 
then has a mean value, of about 0.25, which shows that 
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Figure 6: (a-d): (a-c): Evolution of an orbit as the flattening parameter a changes with time, following the first of equations (21). The orbit 
starts as chaotic and ends as a regular box orbit, in the spherical galaxy, (d-down right): Evolution of the L.C.E of this orbit. Details are 
given in the text. 



the motion has become chaotic. 

The above analysis shows, that the role of the exter- 
nal perturbation is not only to affect the nature of orbits 
(regular or chaotic), but also to change the profile of the 
velocities of the stars. All numerical calculations suggest, 
that as the external perturbation increases and the orbit 
tends to be chaotic, the velocity increases and its profile 
changes. 

5. Discussion and conclusions 

In this work we have used a simple potential, in or- 
der to study the dynamical behavior of a galaxy with an 
additional external perturbation, caused by a companion 
galaxy. Our aim was to investigate the consequences of 
the external perturbation on the character of motion (reg- 
ular or chaotic). Moreover, we have tried to connect the 
strength of the external perturbation, the flattening pa- 
rameter and the radius of the nucleus, with the conserved 
component of the angular momentum L z . 

The numerical calculations suggest, that strong exter- 
nal perturbations cause large chaotic regions on the phase 
plane. Spherical galaxies display smaller chaotic regions 
than flat ones. On the other hand, numerical and the- 
oretical outcomes indicate that linear relationships exist 
between the critical value of the angular momentum L zc 
and the dynamical parameters of the system, that is the 



strength of the external perturbation A, the flattening pa- 
rameter a and the radius of the nucleus c. Moreover the 
extent of the chaotic regions observed in the r — p r phase 
plane, increases linearly, as the strength of the external 
perturbation and the flattening parameter increases. On 
the contrary, the percentage covered by chaotic orbits in 
the phase plane, decreases linearly as the scale length of 
the nucleus increases and becoming less dense. 

The magnitude of the core radius c together with the 
value of the angular momentum L z , are two basic parame- 
ters for the dynamical system to display regular or chaotic 
motion. It was found numerically, that a linear relation- 
ship exists between the critical value of the angular mo- 
mentum and the corresponding radius of the central con- 
centration. The above mentioned relationship can also be 
obtained, using semi-theoretical methods. 

An important role on the evolution of the chaotic mo- 
tion is played by the flattening parameter a of the dy- 
namical system. For a given value of the radius c and the 
external perturbation A, there is a linear relationship be- 
tween a and the critical value of the angular momentum 
L zc . This strongly suggests, that low angular momentum 
stars, display chaotic motion in highly flattened or per- 
turbed elliptical galaxies, having a dense nucleus. Here we 
must note, that this behavior is similar to that observed 
in disk galaxies, studied by Caranicolas & Innanen [3]. 

It is well known, from earlier work, that galaxies with 
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Figure 7: (a-d): (a-b): Evolution of an orbit as the external perturbation A changes with time, following the second of equations (21). The 
orbit starts as regular orbit and gradually becomes chaotic, (c-down left): Evolution of the total velocity of the orbit and (d-down right): 
Evolution of the L.C.E of this orbit. Details are given in the text. 



dense and massive nuclei, produce chaotic regions [4] . The 
phenomenon is strongly connected with the critical value 
of angular momentum L zc and the mass of the nucleus 
M n . The corresponding relationship [L zc ,M n ] was found 
to be linear [3]. Thus, we observe that the behavior of 
low angular momentum orbits in galaxies with massive 
and dense nuclei, is similar to the same orbits in galax- 
ies with strong external perturbation. In both cases we 
observe large chaotic regions on the phase plane and the 
corresponding relationships [L ZCl M n ], [L zc , A], [L ZCl a] and 
[L ZCl c] are linear. 

It is important to note that, in order to observe chaos, 
the external perturbation must be combined with a dense 
nucleus. Numerical calculations not given here suggest 
that, when: A = 15, a = 1.9, L z = 10, no chaotic motion is 
observed when 0.5 < c < 1. Furthermore, it is well known 
that when A = [15], the system displays small chaotic 
regions, only for small values of c. The importance of this 
fact, can be shown as follows. 

The density corresponding to potential (1) is 

,(M = ^ ( ° + 2)c2 + tt < r2 - ( V 2) ^. (22) 
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For small values of r and z, p(r, z) tends to the value 
v 2 (a + 2) 
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Thus, we see that the density near the center increases as 
1/c 2 . Equation (23) shows the critical role of the radius of 
nucleus in the character of motion. Therefore, we can say 
that our numerical experiments suggest that large external 
perturbations can cause significant chaos, if combined with 
high density objects in the central regions of galaxies. 

An estimation of the total mass of the primary galaxy, 
can be obtained, if we assume a spherical galaxy of radius 
R. In this case a = 1 and the density (22) becomes 
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where w 2 = r 2 + z 2 . The mass of the galaxy is 



M G = 4-7T 
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Taking R = 20kpc, c = 0.25kpc we obtain the value: Mq ~ 
4500 mass units, that is 1.05 x 1O 11 M . Furthermore, we 
assume that the two bodies (the primary and the com- 
panion galaxy), are moving in circular orbits around the 
center of mass of the system, with an average period T = 
5 x 10 9 yr. The distance between the primary galaxy and 
its companion is A = 70kpc. The mass of the companion 
galaxy, can be obtained from the Kepler's third law 



Av 2 = G (M G + M c ) , 



(26) 



where v is the relative velocity of the two galaxies. From 
relation (26) and for the given values of the involving pa- 
rameters, we have that M c ~ 670 mass units, which is 
equal to 1.55 x 10 10 M Q . Our numerical results, can now 
be compared with observational data from the binary stel- 
lar system, consists of the giant elliptical galaxy NGC 
3379(M105) and its companion NGC 3384. 

Interesting results are found, in the case where the po- 
tential is time dependent. In this case, the numerical cal- 
culations indicate that the motion evolves from regular 
to chaotic, or from chaotic to regular depending on the 
particular values of the parameters a and A. The main 
conclusion is that, as the external perturbation increases, 
the orbits become chaotic and the total velocity increases. 
Always we must remember, that all numerical calculation 
suggest that a large external perturbation, as described 
by our dynamical model (2), is responsible for producing 
a large amount of chaos, only if the mass density near the 
central regions is high. 

Forty years ago, galactic activity and interactions be- 
tween galaxies were viewed as unusual and rare. Nowa- 
days, they seem to be segments in the life of many galaxies. 
From the astrophysical point of view, in the present work, 
we have tried to connect galactic activity and galactic in- 
teractions with the nature of orbits (regular or chaotic) 
and also with the behavior of the velocities of stars in the 
primary galaxy. We consider the outcomes of the present 
research , to be an initial effort, in order to explore matters 
in more detail. As results are positive, further investiga- 
tion will be initiated to study all the available phase space, 
including orbital eccentricity of the companion and its in- 
clinations to the primary galaxy. 
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